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We show that the Fano representation leads to a particularly simple and appealing form of the 
quantum process tomography matrix \ F , in that the matrix \ F * s real > the number of matrix ele- 
ments is exactly equal to the number of free parameters required for the complete characterization of 
a quantum operation, and these matrix elements are directly related to evolution of the expectation 
values of the system's polarization measurements. These facts are illustrated in the examples of 
one- and two-qubit quantum noise channels. 
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I. INTRODUCTION 

The characterization of physical, generally noisy pro- 
cesses in open quantum systems is a key issue in quan- 
tum information science [1, 2]. Quantum process to- 
mography (QPT) provides, in principle, full information 
on the dynamics of a quantum system and can be used 
to improve the design and control of quantum hard- 
ware. Several QPT methods have been developed, in- 
cluding the standard QPT [2, 3, 4, 5], ancilla-assisted 
QPT [6, 7, 8], and direct characterization of quantum 
dynamics [9] . In recent years QPT has been experimen- 
tally demonstrated with up to three-qubit systems in a 
variety of different implementations, including quantum 
optics [8, 10, 11, 12, 13, 14, 15, 16], nuclear magnetic 
resonance quantum processors [17, 18, 19], atoms in op- 
tical lattices [20], trapped ions [21, 22], and solid-state 
qubits [23, 24]. 

Any quantum state p can be expressed in the Fano 
form [25, 26, 27] (also known as Bloch representation). 
Since the density operator p is Hermitian, the parameters 
of the expansion over the Fano basis are real. Further- 
more, due to the linearity of quantum mechanics, any 
quantum operation p — > p' = £{p) is represented, in the 
Fano basis, by an affine map. 

In this paper, we point out that in standard QPT it 
is convenient to compute the QPT matrix in the Fano 
basis. Such process matrix, x F , has the following advan- 
tages: (i) the matrix elements of \f are re& l an d (ii) the 
number of matrix elements in \ F is exactly equal to the 
number of free parameters needed in order to determine a 
generic quantum operation. Furthermore, the x F -matrix 
elements are directly related to the modification, induced 
by the quantum operation £ , of the expectation values 
of the system's polarization measurements. We will illus- 
trate our results in the examples of one- and two-qubit 



quantum noise. In particular, we will determine in the 
X F -matrix the specific patterns of various quantum noise 
processes. Finally, we will discuss the number of free pa- 
rameters physically relevant in determining a quantum 
operation for a two-qubit system exposed to weak local 
noise. 



II. FANO RAPRESENTATION OF THE 
STANDARD QPT 

To simplify writing, we discuss the Fano representation 
of the standard QPT only for qubits, even though the ob- 
tained results can be readily extended to qudit systems. 
Any n-qubit state p can be written in the Fano form as 
follows [25, 26, 27]: 

P = Jj X! C ai ... an a ai <g) • • • (g) <7 a „, (1) 

Qi , . . .,cx n —x,y,z,I 

where N = 2™, cr x , a y , and <r z are the Pauli matrices, 
oj = 1, and 

c ai ... an = Tr(o- ai ® ••• ®a an p). (2) 

Note that the normalization condition Tr(p) = 1 implies 
c/.../ = 1. Moreover, the generalized Bloch vector b = 
\Pa}a=i,...,N 2 -i is rea l due to the hermiticity of p. Here 
b a = c ai ,.. an , with a = J^fc=i ik^ n ~ k , where we have 
defined ik — 1, 2, 3, 4 in correspondence to — x, y, z, I. 
Note that from 1 to n qubits run from the most significant 
to the least significant. For instance, for two qubits (n = 
2), the N 2 — 1 = 15 components of vector b are ordered 
as follows: 

t* (^1) ^2; •••) ^15) (.Cxxi Cxyi Cxzi Cxi: Cyxi ^yyj 

Cyz 3 7 , C zx , C Z y , C zz , C Z J , Cj x , Cjy , Cj z ) . 
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Due to the linearity of quantum mechanics any quan- 
tum operation p — > p' = £(p) is represented in the Fano 
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basis {<J ai 



<7 an } by an affinc map: 
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(4) 



where M is a (A 2 — 1) x (TV 2 — 1) matrix, a a column 
vector of dimension A 2 — 1 and the null vector of the 
same dimension. 

All information about the quantum operation £ is con- 
tained in the A 4 — A 2 free elements of matrix A4, namely 
in the matrix 



X F = M 



(5) 



To obtain the QPT matrix x F from experimental data, 
one needs to prepare A 2 linearly independent initial 
states {pi}, let them evolve according to the quantum 
operation £ and then measure the resulting states {p- = 
£(pi)}. If we call TZ the N 2 x A 2 matrix whose columns 
are given by the Fano representation of states pi and TZ' 
the corresponding matrix constructed from states p' i7 we 
have 



and therefore 



TZ' = MTZ, 



m = n"R- 



(6) 



(7) 



As it is well known [2], the standard QPT can be per- 
formed with initial states being product states and local 
measurements of the final states. As initial states {pi} we 
choose the 4™ tensor-product states of the 4 single-qubit 
states 



|0>, 



|1>, ^=(|0> + |1», -^(|0)+i|l)). (8) 

To estimate TZ' , one needs to prepare many copies of each 
initial state pi, let them evolve according to the quantum 
operation £ and then measure observables cr ai <£>• • -®<J an . 
Of course, such measurements can be performed on the 
computational basis {|0), provided each measure- 

ment is preceded by suitable single-qubit rotations. 



III. SINGLE-QUBIT SYSTEMS 



The matrix TZ corresponding to basis (8) reads 



TZ 



10 
1 
1-10 
1111 



(9) 



Therefore, 



n- l = 



_i _i l i 

2 2 2 2 

_1 _1 _1 1 

2 2 2 2 

10 











(10) 



The coefficients (c x ,c yi c z ) in the Fano form (1) are 
the Bloch-vector coordinates of the density matrix p in 
the Bloch-ball representation of single-qubit states. We 
need A 4 — A 2 = 12 parameters to characterize a generic 
quantum operation acting on a single qubit. Each pa- 
rameter describes a particular noise channel (like bit flip, 
phase flip, amplitude damping,...) and can be most con- 
veniently visualized as associated with rotations, defor- 
mations and displacements of the Bloch ball [1, 2, 28]. 
Here we point out that these noise channels lead to spe- 
cific patters in the state process matrix x F ■ 

For instance, for the phase-flip channel, 

p' = £{p)=pa z pa z + (l-p)p, (0< P <i), (11) 



we have 



TZ' 



l-2p 

1 - 2p 
1-10 
111 1 



(12) 



We can then compute M = TZ' R 1 , and the first three 
lines of M correspond to the state matrix 



y(pf) 



l-2p 
l-2p 
10 



(13) 



Therefore, the Bloch ball is mapped into an ellipsoid with 
z as symmetry axis: 



(1 - 2p)c x 

(1 - 2 P )Cy 



(14) 



As a further example, we consider the amplitude damp- 
ing channel: 



(15) 



fe=0 



with the Kraus operators 



i?o = |0)(0| + v ^p|l)(l|, E 1 = y /p\0)(l\, (0<p<l). 

(16) 

In this case we obtain 



,(ad) 



VT^p 
V^P 
1 -p p 



(17) 



The Bloch ball is deformed into an ellipsoid, with its 
center displaced along the z-axis: 



c 'x = V 1 -P c x, 

c 'y = ^J^-Pcy, 

c 'z = (! -P)c z +P- 



(18) 
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IV. TWO-QUBIT SYSTEMS 



Matrices TZ and 1Z 1 corresponding to the 16 tensor-product states of singlc-qubit states (8) read as follows: 
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(19) 
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The coordinates {c aia2 } in the Fano form (1) are 
the expectation values of the polarization measurements 
{u ai ® &a 2 }- The coefficients in the state matrix x F 
representing a quantum operation £ can therefore be in- 
terpreted in terms of modification of these expectation 
values. 

For instance, let us assume that the two qubits are in- 
dependently exposed to pure dephasing, that is, to quan- 
tum noise described by the phase- flip channel (11), with 



the same noise strength p for both qubits. The process 
matrix for such uncorrelated dephasing channel is given 
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by 



cess matrix 



Y (ud) = 

A p 
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where g = 1 — 2p. Correspondingly, the mapping for 
the expectation values of the polarization measurements 
reads 



— y °o 



(22) 



where rrii = 1 for on = x,y and rrii — for on = z, I. 

As an example of nonlocal quantum noise, we consider 
a model of fully correlated pure dephasing. We model 
the interaction of the two qubits with the environment 
as a phase-kick rotating both qubits through the same 
angle 9 about the z axis of the Bloch ball. This rotation 
is described in the {|0), |1)} basis by the unitary matrix 



R z (0) 



e-^ 




e - i f 


e l i 




e 4 i 



(23) 



We assume that the rotation angle is drawn from the 
random distribution 



Y (cd) = 

A F 



h 00 fcOOOOOOOOO 

h 00 -k 0000000000 

g 0000000000 

g 0000000000 

0-fc00 h 0000000000 

k 00 ZiOOOOOOOOO 

00 9 

o o oo o oo^ooooooo 

00 OOOjOOOOOO 
00 j 
00 0000010000 
00 0000001000 
00 OOOOOOOjOO 
00 00000000 5 
00 0000000001 

(26) 

i=l(l+ 5 4 ), k = \(1 - g 4 ). It is clear 
that the process matrix (26) for correlated dephasing has 
a pattern that allows to clearly distinguish it from the 
process matrix (21) for the uncorrelated dephasing. 

It is also obvious that, if there exists partial previous 
knowledge of the dominant noise sources, it is not nec- 
essary to construct the whole state process matrix \ F 
in order to characterize the quantum operation. For in- 
stance, if we know a priori that dephasing is the main 
source of noise and we wish to estimate its degree of 
correlation, it is sufficient to prepare, for instance, the 
initial state p — |(|0) + |1))® 2 and measure the x- and 
y-polarizations of both qubits for the final state p' . The 
initial state is fully polarized along x, and therefore 



(21) where g = e 



Cxx 1? 

c yy = 0- 



(27) 



For the final state, in the case of fully correlated dephas- 
ing 




)(cd) =h=\[l 
)^)=k = \[l 



9 



(28) 



p(9) 



-e 4A . 



(24) 



while the expectation values of the xx- and yy- 
polarization measurements are remarkably different for 
uncorrelated dephasing: 



Therefore, the final state p', obtained after averaging over 
6, is given by 



/+OC 
de P (6)R z (e)pRl(0). (25) 
-00 



For this correlated dephasing channel we obtain the pro- 



(c' x J^=g 2 , 

KJ (ud) = 0. 



(29) 



While in general two-qubit quantum operations de- 
pend on iV 4 — N 2 — 240 real parameters, an important 
question is how many parameters are physically signifi- 
cant. The answer of course depends on the specific noise 
processes. However, a clear answer can be given assuming 
that external noise is weak and local, that is to say, it acts 
independently on the two qubits. In this case, local noise 
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is described by 24 parameters, 12 for each qubit. Undc- 
sired coupling effects (cross talk) between qubits can be 
characterized with only three additional parameters, 6 X , 
9 y , and 8 Z . Indeed, any two-qubit unitary transformation 
U can be decomposed as [29, 30, 31] 

U = {A 1 ® B 1 )e i ( e *' 7x ®' 7 * +e v a v®' 7 v +9 *' 7 *®' 7z \A2 <g> B 2 ), 

(30) 

with Ai, A 2 , B\, and B 2 appropriate singlc-qubit uni- 
taries. In the limit of weak noise, the state matrix x F 1S 
simply given by the sum of the contributions of each noise 
channel. Therefore, the local unitaries A\, A 2 , B\, and 
B 2 only change the 24 local noise parameters and overall 
we need 24 + 3 = 27 <S 240 parameters to describe the 
quantum noise. In the symmetric case in which the local 
noise parameters are the same for both qubits the number 
of free parameters further reduces to 12 + 3 = 15. The 
above argument can be easily extended to many-qubit 
systems. Due to the two-body nature of interactions, we 
need to determine TV = 12n+3 "^" 2 ~ 1 ^ parameters to char- 
acterize noise. Note that TV = O(logTV) «JV 4 - TV 2 . Of 
course, cases with strong or nonlocal noise would require 
a larger number of free parameters. 



V. CONCLUSIONS 



We have shown that the Fano representation of the 
standard QPT is convenient, since the process matrix % F 
is real and the number of matrix elements is exactly equal 
to the number of free parameters required for the com- 
plete characterization of a generic quantum operation. 
Moreover, the matrix elements of x F are directly related 
to the evolution, induced by the quantum operation, of 
the system's polarization measurements. We have also 
shown that quantum noise channels have specific patterns 
in the Fano representation of \ F . Finally, we have shown 
that in the case, of interest for quantum information pro- 
cessing, of weak and local noise the number of relevant 
noise parameters is TV = O (log TV) < TV 4 - TV 2 , that 
is, much smaller than the number of parameters needed 
to determine a generic quantum operation. In this case, 
the x F -matrix is very sparse and therefore the number 
of polarization measurements needed to reconstruct it is 
much smaller than for a generic quantum operation, thus 
considerably reducing the QPT complexity. 
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